Nonclassical effects and off-diagonal couplings in a model for FeBr 2 
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Using Monte Carlo techniques, we show that the re- 
cently experimentally observed transitionlike phenomena in 
the transverse spin-ordering close to the anomalies in the an- 
tiferromagnetic phase of FeBr2 (Petracic et al. [Phys. Rev. 
B 57, R11051 (1998)]) may result from the quantum nature 
of the S=l spins and an off-diagonal exchange between axial 
and planar spin components. 



Recently, Petracic et al.0 presented experimental ev- 
idence for a weakly first-order phase transition in the 
antifcrromagnetic phase of FeBr2, when applying a mag- 
netic field under non-vanishing tilting angle with respect 
to the c-axis of the hexagonal crystal, extending previous 
work with an axial field.BH They argue that the measured 
jumps in the magnetization parallel and perpendicular to 
the field are related to the surprisingly sharp peaks near 
the anomalies in the specific heat which had been mea- 
sured before in axial magnetic fields.Q 

They suggest that quantum effects in a S=l anisotropic 
Heisenberg magnet and off-diagonal exchange between 
axial and .planar spin components allowed by the crystal 
symmetryd may be crucial in explaining the transitionlike 
phenomena. 

Indeed, our analyses on an appropriate model support 
their suggestion. In particular, we considered the follow- 
ing S = 1 Hamiltopian, being believed to describe quite 
realistically FeBrJjllJ, 
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The first term describes anisotropic, a < 1, exchange in- 
teractions, Jy, between nearest (Ji), next-nearest (J2), 
and third (J3) neighbor spins in the triangular layers per- 
pendicular to the c-axis and between spins in adjacent 
layers connected by equivalent exchange paths (J') .ETcHl 
The second term describes the single-ion anisotropy of 
strength D. The third term refers to the tilted field, hav- 
ing components both perpendicular (H z ) und parallel 
(Ha,) to the triangular planes. The last term H d ( Jxz) 
describes the off-diagonal bilinear couplings between ax- 
ial (S z ) and planar (S x and S y ) components of neigh- 
boring spins in the layers, with the strength J^-lI For 
instance, for neighboring spins with the same value of y, 
one has 

H'od = —Jxz (S-Sj + S^Sj) . (2) 

The corresponding terms for other pairs of nearest- 
neighbor spins follow from appropriate translations and 
rotations.cl 



We performed Monte Carlo simulations on that Hamil- 
tonian in the quasi-classical approximation, i.e. fixing 
the spin length to be y/2, with S z = ±1,0. The exchange 
constants Jij and the anisotropy factoiLo; were chosen as 
beforeEHl (based on spin-wave analysestm) , with the ten- 
fold coordinated antifcrromagnetic couplings J' between 
adjacent triangular layers, with nearest-neighbor, Ji= - 
16.75 J', and competing third neighbor, J 3 = -0.29 J1, 
interactions in the planes perpendicular to the c-axis, 
and with an Ising-type anisotropy factor a = 0.78. We 
then varied the degree of the Ising-like anisotropy due 
to the single ion term, D, the axial, H z , and planar, IP, 
components of the external magnetic field, as well as the 
strength, J XZl of the off-diagonal exchange between the 
axial and planar spin components. 

In an axial field, H x = and H z being close to the crit- 
ical field, H z o = —2 J' (setting the Boltzmann constant 
equal to one), at zero temperature, the peculiar shape 
of the temperature dependent specific heat C(T) found 
experimentallyu'EJ can, indeed, be reproduced by tuning 
the remaining parameters of the model, as shown in Fig. 
1. In particular, the maximum at the lower tempera- 
ture shows a pronounced shoulder, associated with the 
anomaly in the z-components of the spinsm-I originating 
from the relatively weak, albeit dominating, intralayer 
couplings as compared to the interlayer couplings_as well 
as from the high interlayer coordination numberEI til. The 
superimposed peak shifts towards the critical tempera- 
ture, T c , and gets sharper as 3 XZ increases. That peak, at 
T X y, reflects the disordering of the antifcrromagnetic low 
temperature state of the xy- or planar components of the 
spins, being aligned ferromagnetically in each in the tri- 
angular planes, with an antiferromagnetic arrangement 
between subsequent layers. The corresponding disorder- 
ing of the z- or axial components occurs at T c . The two 
processes are largely decoupled due to the quantization 
of the spins. Note that the anomaly and the transition- 
like peak are not intimately related, and their positions 
can be moved relatively to each other by changing, for 
instance, the off-diagonal exchange 3 XZ and the degree 
of Ising-like anisotropy D. 

Applying a non-axial field, H x > and keeping H z 
fixed, the xy-components of the spins order in the an- 
tiferromagnetic phase, at low temperatures, in a spin- 
flop state, in which the x-components of the magnetiza- 
tion per layer are almost identical in all layers, but the 
y-components of the magnetization in adjacent planes 
have opposite signs. Increasing the temperature, the 
y-component of the magnetization per layer changes 
rapidly at the temperature T xy , well below the transi- 
tion temperature T c of the antiferromagnetic phase, see 
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the Monte Carlo data depicted in Fig. 2. The change 
in the y-component may be associated with a first-order 
transition, characterized by a jump in the magnetiza- 
tion per layer, when the off-diagonal exchange 3 XZ ex- 
ceeds a critical value (at about J xz = 18 J'). Other- 
wise, one finds a-drastic, but presumably analytic change. 
Experimentally!!], one is at the border of these two scenar- 
ios, observing a 'weakly first-ordecl transition at T xy — 
Note that the experimental findingsEI are compatible with 
the spin-flop state. 

In summary, our simulational data show that inclu- 
sion of the off-diagonal couplings in the spin components 
allowed by the crystal symmetry and of the quantum 
nature of the S=l spins in the Hamiltonian for FeBr2 
is, indeed, sufficient to describe the transitionlike fea- 
tures in the specific heat and the magnetization observed 
experimentally recently. The features are related to a 
rapid change or transition in the planar spin components, 
which, in turn, seems to occur rather fortuitously closely 
to the anomaly of the axial spin components in the an- 
tiferromagnetic phase of FeBr2 discussed extensively be- 
fore. 
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FIG. 1. Temperature dependent specific heat C as obtained 
from Monte Carlo simulations of Hamiltonian (1) with J X z = 
16.2 J', D= -8.1 J', and H z = -1.8 J' (= 0.9H z0 ), for systems 
of 20 3 spins. 
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FIG. 2. Temperature dependence of the x-, y- and 
z-components of the magnetization per layer in odd (full sym- 
bols) and even (open symbols) planes, using the same param- 
eters as in Fig. 1, with an additional planar component of the 
field, H x — 0.75 H z , simulating systems of 30 3 spins. 
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